Fermi surface topology and ferromagnetic superconductivity in UGe2 
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We consider a Stoner ferromagnet in presence of a quasi-one dimensional Fermi surface in its spin 
majority band. Assuming a twin 5-function peaked density of state due to the low dimensionality, we 
computed the single particle self-energy. There appears an additional divergence in the self-energy 
and hence in the effective mass inside the ferromagnetic phase (besides the standard logarithmic 
divergence at the Stoner critical point). Since such an effect is purely due to density of states, it 
might correspond to a first order phase transition making the ferromagnetic phase into two distinct 
phases. This result is in qualitative agreement with the recent specific heat capacity measurement. 
We also discuss its relevance to the superconducting state in UGe2. 

PACS Numbers: 71.27.+&, 74.20.Mn, 75.10.Lp, 75.50. Cc 
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Superconductivity (SC) and ferromagnetism (FM) are 
two very different cooperative phenomena and the ques- 
tion regarding their coexistence is very important theo- 
retically as well as experimentally. Due to their antago- 
nistic character, it was believed that they mutually act 
against each other. But the recent discovery of FMSC 
in UGe2 [1], URhGe [2] and ZrZna [3] has ruled out this 
possibility. The behaviour of these materials is an exam- 
ple of a more general phenomenon where a novel state 
appears on the boarder of magnetism at low tempera- 
ture. In what follows, we concentrate on the material 
UGe2, where SC is entirely covered within the FM phase 
and disappears in the paramagnetic (FM) region. 

UGe2 is an itinerant FM [4] with a Curie tempera- 
ture of Tc= 52 K at ambient pressure. Upon increasing 
pressure, it becomes SC with a maximum transition tem- 
perature of Ts ^ 0.8 K in the pressure range of 1 GPa 
to 1.7 GPa. Furthermore, there appears an additional 
phase line T^; [4] inside the FM phase (see Fig.l) divid- 
ing it into two distinct phases (FMl and FM2). This has 
been inferred from the strong anomaly seen in the resis- 
tivity measurements [4] , a small enhancement in the spe- 
cific heat capacity [5] , lattice expansion [6] and a change 
in the character of the Fermi surface measured in the de 
Haas van Alphen experiments [7] . This might be thought 
to be associated with the formation of charge and spin 
density waves (CSDW) since the band structure calcula- 
tions [8,9] indicate that this material is prone to nesting. 
Until now, no density fluctuations has been observed in 
the neutron scattering measurements. 

It is already known that in an itinerant weak FM, the 
magnetic fluctuations are enhanced in the vicinity of the 
quantum phase transition and it gives rise to triplet SC 
[10]. The standard phase diagram proposed in such a 
formalism is the appearance of SC dome both in the PM 
as well as FM phases close to the Curie transition. The 
SC transition temperature estimated in such a formula- 
tion is comparable in both the phases. Recently, similar 
results have been obtained by solving the Eliasberg equa- 
tions rigorously [11]. However, it has also been pointed 



out that SC in the FM region can be enhanced at least 
50 times than that of the PM phase [12]. This can be 
understood due to the mode coupling among transverse 
magnons and longitudinal spin fluctuations in the FM 
phase which is absent in the PM region. However, the 
FM state of UGe2 is so magnetically anisotropic that the 
presence of magnons seems an unlikely explanation for 
the enhancement of SC in the system. Furthermore, in- 
spired from the band structure calculations [8,9], Watan- 
abe and Miyake [13] postulated the existence of CSDW 
inside the FM phase, in analogy with the a-phase of ura- 
nium. In such a formalism, CSDW fluctuations at high 
wave vector couple to the magnetization in such a way 
that it yields a mass enhancement near CSDW transi- 
tion. It also provides a pairing mechanism for SC in this 
system. 
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FIG. 1. Schematic T-P phase diagram in UGe2 showing 
the phase line which separates the two magnetic phases 
FMl and FM2. 



Very recently, the pressure dependence of magnetiza- 
tion measurement at low temperature by Pfleiderer et 
al., [14] indicates that the two transitions (Curie as well 
as the transition between FMl and FM2) are of first or- 



1 



dcr in nature. They further pointed out that even if 
the low temperature uniform longitudinal susceptibility 
undergoes a large change between FMl and FM2, it be- 
comes almost pressure independent within each phase. 
This suggests that the SC is not driven due to its prox- 
imity to a FM quantum critical point; rather, it could 
be associated with a sharp spike in the electronic den- 
sity of states (DOS). This could also be responsible for 
the low pressure magnetic transition. Keeping this in 
mind, Sandeman et al. [15] studied the zero tempera- 
ture Stoner model in presence of a two Lorentzian peaked 
DOS phenomenologically. They minimised the total en- 
ergy density of the system with respect to the magneti- 
zation and obtained the magnetization-Stoner exchange 
(M-U) phase diagram. For suitable choices of the DOS 
parameters, they showed that one can obtain two first 
order transitions. They also discussed its relevance to 
superconductivity. 

In this letter, we study the effect of changing Fermi 
surface topology in an itinerant FM. Due to the hidden 
quasi-one dimensionality, we assume a double 5-function 
peaked DOS and compute the single particle self-energy. 
Besides the standard logarithmic divergence due to spin 
fluctuations at the Stoner critical point, it yields an an- 
other divergence inside the FM phase dividing it into two 
distinct phases. Since this is purely due to DOS effect, 
one might argue this phase line to be of first order in na- 
ture. The effective mass at this transition gets enhanced, 
which is in qualitative agreement with the recent specific 
heat capacity measurement. We also discuss the role of 
such DOS in the superconducting phase of UGe2. 

Since the present work is motivated from the material 
UGe2 where the 5f electrons from the U-atoms play dual 
roles both for FM and SC, we can start with a model 
which has Hubbard type exchange interaction. This is 
inferred from the specific heat data [4] which is at least 
ten times lower than the standard heavy Fermion ma- 
terials. Further, the measured low residual resistivity 
(0.2 /ifi) yields a carrier mean free paths of several thou- 
sand angstroms [1]. Thus, the 5f electrons from the 
U-atoms arc considered to be itinerant but of course, 
strongly correlated. We presume that the relevant mag- 
netic behaviour of the system is described by Stoner 
RPA-mean field theory [16]. Further, due to the pro- 
nounced magnetic anisotropy of UGe2, and as a conse- 
quence, the absence of tranverse spin modes makes the 
Stoner approach better at low temperature. The energy 
for single particle excitation in a Stoner FM is given 
as, ek,(T = €k — H + aAp. Here, cr = ± for spin up 
(down) bands, Ap is the order parameter (imiform mag- 
netization) which characterizes the itinerant FM state by 
2 Air = J7(<nj^ > — < n-i >), U is the on-site Hubbard 
interaction energy and nk,cr = c[ ^Ck,a is the density op- 
erator with wave vector k and spin projection a. The 
occupation of each spin band a is, Tia = '°(^)'^^' 



where {—W/2) is the bottom of the band and spin a oc- 
cupies energy states up to /Uo-- The DOS here is given 
by p(e). We also consider the total number of spins, 
N = n^' + to be fixed. Thus, the chemical potential 
/Lto- of each spin band is completely determined by the 
particular magnetization and the electron number. 

The single particle propagator for a weak FM is given 
as, 
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where, the residue Z, contains informations about the 
quasiparticle properties (which has to be less than one for 
well-defined quasiparticles as in the framework of Fermi 
liquid theory). The incoherent part of the Green func- 
tion (Ginc), which describes the high energy processes, is 
a smooth function of iv and k. In such a situation it can 
only renormalize the properties at the Fermi surface with- 
out introducing any new physics. Such a liquid has two 
types of low energy collective spin excitations: transverse 
magnons and longitudinal spin fluctuations. For UGe2, 
due to the large magnetic anisotropy, we consider only 
the latter in what follows. The relevant expression for 
the longitudinal spin fluctuation [17,18] is given as, 
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Here, p{0) is the DOS at the Fermi level in the FM state, 
ry = (A — 1) ~ A = Up{0), A and C are the param- 
eters related to ep and p(0). Now, using the standard 
many body technique [19], one can calculate the single 
particle self-energy in the spin majority band as, 

A^ 

T,{k,u) = -g X] G'^{k+p,iUn + ipn)xi{P,Wn), (3) 

where Hujn represents Matsubara frequencies. At this 
stage of the calculation, we would like to point out that 
the wave vector sum in the above equation can be per- 
formed by using the DOS of the FM state. The standard 
way to get a proper form of the DOS is to analyze the 
band structure data. We do exactly this for the present 
case. It has already been noted from the band struc- 
ture calculations [8,9] that the Fermi surface in UGe2 
is quasi-cylindrical. In presence of exchange splitting, 
the minority electron surface contracts while the major- 
ity spin counterpart expands towards the boundary of 
the Brillouin zone and gets cut-off by the zone wall on 
both sides. This produces two large and roughly parallel 
sheets reminiscent of a quasi-one dimensional system. It 
has been known from text books that in a perfectly one- 
dimensional tight binding system, the DOS has square 
root singularity and it diverges at both the band edges. 
One can show that by including the higher order har- 
monics in the dispersion, it is possible to get double peak 
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structure in the DOS as has been shown by Sandeman et 
al [15]. 

To make the p-sum tractable analytically in equation 

(3), wc assume a following form of DOS, which has dou- 
ble (5-function peaked structure, written as. 



p(e)=p(0)[l+a(5(e-6) + 5(e + d))]. 



(4) 



Here, 'a', '6' and 'd' are the parameters which determine 
the strength and the positions of the DOS. These param- 
eters can be determined from the band structure calcula- 
tions. Moreover, it should be noted that the energy e in 
the DOS is measured with respect to the Fermi energy of 
the spin majority band (fFt)- In what follows, we show 
that this hidden quasi-one dimensional character associ- 
ated with the large electron sheet of the majority spin 
Fermi surface which yields the above DOS (4), will have 
profound consequences. 

Now, by performing the Matsubara frequency summa- 
tion as well as the p-sum by using the above form of the 
DOS, one obtains, 
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rj + Ab 



7]- Ad 
(5) 



where P = Akp and Q = a/ep are the dimensionless 
variables. Here, A = P/ep = Akp/ep. By differentiat- 
ing the real part of the self-energy, one can compute the 
effective mass, which turns out to be, 



parameter (ry), we don't go for a quantitative fit because 
the experiments are done with respect to pressure. We 
would like to further make some comments on the second 
term in equation (6). From the pressure dependence of 
magnetization, it is clear that even the Stoner transition 
is first order. In a standard weak Stoner FM, the PM-FM 
transition is second order. Thus, we can make a choice, 
6 = in the equation (6) so that both the first and the 
second term in the above equation cause mass divergence 
at at the same point, 77 = 0. Thus, one of the peak in 
the DOS should lay exactly at the Fermi level of the spin 
majority band (epi) (see Fig.2). 




FIG. 2. Single particle DOS whose one of the peak coin- 
cides exactly at the Fermi level of the spin majority band 
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(6) 



This is the central result in the manuscript. Recalling the 
weak itinerant FM, it is known that the effective mass as 
well as the specific heat capacity diverges logarithmically 
near the Stoner critical point (r; ^ 0) . This is exactly the 
first term in the above expression whereas the next two 
terms are due to the double (5-function peaked DOS. It 
is obvious from the last term in the above equation that 
there is an algebraic divergence for finite 77, which causes 
a phase transition inside the FM phase. This corresponds 
to the FM1-FM2 transition which has been discussed ear- 
lier. The specific heat capacity measurement by Tateiwa 
et al., [5] shows that the value C/T increases approx- 
imately to 95 mJ/mole K"^ at pressure about 1.15GPa 
which is nearly the pressure where FM1-FM2 transition 
occurs. This causes an enhancement of the the effective 
mass as much as three times than that at ambient pres- 
sure. We can obtain such a mass enhancement from the 
above equation by considering suitable parameters. In 
the present formulation, since the enhancement in the ef- 
fective mass is with respect to the dimensionless coupling 



Next, let us consider the effect of such a twin peak DOS 
on SC. Due to the strong magnetic anisotropy in UGe2, 
it is believed that SC in this system is not only equal spin 
pairing (ESP) but also non unitary [20]. Following Fay 
and Appel [10], we can write down the self- consistent gap 
equation as, 



Ar'(fc,ia;) = i^y-'<^'(fc,za;;p,ii/) 



(7) 



where V'^''^ is the irreducible pairing potential. For a 
general /-state ESP case, one can define a BCS type pair- 
ing parameter gf as, gf = p„{0)Vi^''^ , with 
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F-'-(|g|,a; = 0). (8) 



Here, q = fc — p with \k\ « |p| « fc^ and (k-p) = kpCosO, 
q = 2 kp \sin{9/2)\. We can compute the pairing po- 
tential V^'"^ due to longitudinal spin fluctuations. For 
the case of ESP, the ladder diagram will be absent due 
to the Pauli exclusion principle and the contribution to 
V^'"^ will only be given by summing the odd number of 
bubble diagrams. Within RPA, one obtains, 
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1 - u^x'^'x: 



(9) 



where x'^'^ is the particle-hole bubble for a particular spin 
projection a. Computing V°'''^ for a weak itinerant FM 
as has been done earlier [10], one can calculate the BCS 
parameter by doing angular integration. Moreover, in 
the present case, we can compute = gis (we drop the 
cr index since we work in a spin majority band, the sub- 
script 'I'represents the angular momentum I = 1, and 's' 
for the SC state) in presence of a double peaked DOS, 
which turns out to be. 
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where Pg = A/6, Qs = 2a/5eF and As = X/lbep- Now, 
one can solve the self-consistent gap equation within the 
standard BCS approximation which yields, 



Ub Tg = lA4wc exp[-l/fife//], 



(11) 



where geff = gis/{^+gz), is the cfFcctivc coupling and u)c 
is the cut-off due to longitudinal spin fluctuations. It has 
been shown earlier by Brinkman and Engelsberg [21] that 
uic is in fact proportional to rj. Of course, the derivation 
of the SC transition temperature is admittedly not very 
satisfactory due to the fact that the effective pairing in- 
teraction and the self-energy in an weakly FM system are 
energy dependent. So the most elegant method for tack- 
ling this problem is to go for a conventional Eliasberg 
strong coupling analysis. However, the weak coupling 
BCS like formulation provides a qualitative understand- 
ing which is presented above. 

It is clear from equation (6) and (10) that the effec- 
tive mass {gz) as well as the BCS coupling parameter 
{gis) for Z = 1 ESP behaves exactly in the similar way, 
i.e., they are peaked at the Stoncr threshold (j] = 0) and 
at the FM1-FM2 transition (77 ~ Agd). Furthermore, the 
effective coupling parameter (geff) also has similar fea- 
ture as that of gz and gis but the absolute magnitudes 
are different. This indicates that the factor exp[— 1/9,./^] 
also attains a maximum at 77 = and 77 ^ Agd. How- 
ever, due to the presence of the cut-off as a prefactor 
in the expression for fc^ T,, the SC transition tempera- 
ture vanishes at the Stoner threshold while it becomes 
maximum near the FM1-FM2 transition. This is what is 
exactly observed in the material UGe2. Thus, as far as 
the transport and the SC is concerned, the present formu- 
lation qualitatively explains the experimentally observed 
spectrum. 

To conclude, we summarize the main results of the 
present work. We study the role of hidden quasi-one di- 
mensional Fermi surface in the spin majority band of a 
weak Stoner FM. This was inferred from the band struc- 
ture calculations in the material UGe2. To make the 



calculations analytically tractable, we assume a double 
(5-function peaked DOS and computed the single parti- 
cle self-energy. We showed that the effective mass and 
hence the specific heat capacity acquires an additional 
divergence besides its standard logarithmic divergence at 
the Stoner critical point. This ultimately divides the FM 
phase into two distinct phases. Since this is purely due 
to the DOS effect, one might argue this transition to be 
of first order in nature. It is also shown that such a DOS 
yields maximum SC transition temperature at a point at 
which the specific heat capacity gets enhanced. This is in 
qualitative agreement with the phase diagram obtained 
for UGe2. 
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ing the manuscript. 
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